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Abstract. We numerically investigate the ground-state phase diagram of an S = 2 quantum 
spin chain with the XXZ and on-site anisotropics described by H = y\ (S^Sj+i +S^S^ +1 + 
ASjSj +1 ) + D^2j(Sj) 2 , where A denotes the XXZ anisotropy parameter of the nearest- 
neighbor interactions and D the on-site anisotropy parameter. We restrict ourselves to the 
A > and D > case for simplicity. Our main purpose is to obtain the definite conclusion 
whether there exists or not the intermediate- D (ID) phase, which was proposed by Oshikawa 
in 1992 and has been believed to be absent since the DMRG studies in the latter half of 1990's. 
In the phase diagram with A > and D > there appear the XY state, the Haldane state, 
the ID state, the large- D (LD) state and the Neel state. In the analysis of the numerical data 
it is important to distinguish three gapped states; the Haldane state, the ID state and the LD 
state. We give a physical and intuitive explanation for our level spectroscopy method how to 
distinguish these three phases. 



1. Introduction 

In this paper, using mainly numerical methods, we [1] investigate the ground-state phase diagram 
of the S = 2 quantum spin chain described by the Hamiltonian 

H = Y^iSjSj.. + S]S] +1 + AS|5f +1 ) + DJ2(S*) 2 , (1) 

j j 

where A and D are, respectively, the XXZ anisotropy parameter of the nn interactions and the 
on-site anisotropy parameter. We restrict ourselves to the A > and D > case for simplicity. 

The ground-state phase diagram of this model was first discussed by Schulz [2]. In his phase 
diagram for A > and D > 0, there appear the Neel phase, the XY phase, the Haldane phase 
and the large-D (LD) phase. In 1992 Oshikawa [3] predicted, for S > 2 integer quantum spin 
cases, the existence of the intermediate- D (ID) phase, the valence bond picture of which is 
depicted in Fig. 1(b). After that, however, by use of the density-matrix renormalization-group 



Figure 1. Valence bond pictures for (a) the 
Haldane state, (b) the ID state and (c) the LD 
state. Big circles denote 5 = 2 spins and dots 
5 = 1/2 spins. Solid lines represent valence 
bonds (singlet pairs of two 5 = 1/2 spins, 
(1/V2)(t4- - it))- Two 5=1/2 spins in 
rectangles are in the (5t t, 5j* ot ) = (1, 0) state 
and similarly four 5 = 1/2 spins in squares are 
in the (5 to t, 5 t 2 ot ) = (2, 0) state. 

(DMRG) calculation, Schollwock et al. [4, 5] and Aschauer and Schollwock [6] concluded the 
absence of the ID phase. By use of the level spectroscopy (LS) analysis of numerical results of 
exact-diagonalization calculations, Nomura and Kitazawa [7] showed in the case of A = 1 that, 
with the increase of D from zero, the ground state changes as Haldane state =>- XY state =4> 
LD state. The first and second transitions occur at D c \ = 0.043 and D C 2 = 2.39, respectively. 
Since these works it has been believed for a long time that the ID phase does not exist in the 
phase diagram of the present 5 = 2 model. 

The DMRG is a very powerful method for spin chains, especially when the magnitude of the 
spin gap is targeted. However, it is difficult to deal with the phase transition in some cases, 
because the phase transition point is determined by the extrapolated values of the finite-size 
being equal to zero. Thus, it is somewhat hard to obtain accurate results if the magnitude of 
the gap is very small (i.e. the correlation length is very long). This difficulty of zero-or-finite 
problem is not special to the DMRG method but is common to almost all numerical methods. 
On the other hand, the LS method is conspicuous on this point because the phase transition 
point is determined from the crossing point of two related excitations, which is free from the 
zero-or-finite problem. 

The LS method [7, 8] is firmly based on the effective Hamiltonian, renormalization group 
method and conformal field theory. In this paper we give a physical and intuitive explanation 
for our LS method how to distinguish the Haldane, ID and LD states sketched in Fig.l. 

2. A very simple example: 5 = 1/2 chain with bond alternation 

Before discussing 5 = 2 chain problem, let us visit a very simple example of an 5 = 1/2 chain 
with bond alternation described by 

UU8) = Et 1 + (-!) j+1 <^ • S j+1 , 5 = 1/2 (2) 

3 

where 5 is the bond-alternation parameter (— 1 < 8 < 1). For 5 = the ground state of the 
Hamiltonian (2) is the Tomonaga-Luttinger liquid state. We treat an N = 8 system for a while. 
When 5 = ±1 the ground states of the Hamiltonian (2) under the periodic boundary condition 
(PBC) are trivial 

< BC (5 = +1) = [1,2][3,4][5,6][7,8] (3) 
< BC (5 = -1) = [2,3][4,5][6,7][8,1] (4) 

where = (l/v / 2)(tilj — Utj) = — Lm]. It is well known that the ground state for 5 > 0, 
^ BC (5 > 0), is similar to ^ BC (5 = + 1 ) and < BC (^ < 0) is similar to ^j pbc (5 = -1). The 
quantum phase transition occurs at 5 = 0. The state i()q BC (5 = ±1) is the highest energy 
states at 5 = respectively. However, since the crossing of the ground state energy does 




not occur when 8 is swept, as depicted in Fig. 2, we cannot determine the transition point from 
the crossing. Because V PBC (<5 = ±1) are not distinguished by the eigenvalue of the discrete 
symmetry, there occurs the mixing of these two states resulting in the level repulsion. In other 
words, the degeneracy is lifted by the "perturbation". 
Let us impose the twisted boundary condition (TBC) 



S 8 ■ Si =► -Sf Sf - S%Sf + SISf . 



(5) 



In this case the ground state of Hamiltonian (2) for 5 = ±1 cases are, respectively 



V iBO (5 = +l) = [l,2] [3,4] [5,6] [7,8] 
^ TBC (5 = -1) = [2,3] [4,5] [6, 7] {8,1} 



(6) 
(7) 



where {i,i} = (1/ \/2) (Ti4-i + Ijti) = {j^}- Here we introduce the lattice inversion operator X 
about the bond center of the boundary, which works 1 ^ 8, 2 H 7, 3 f* 6, 4 f> 5. From 



IVo TBC (<5 = +l) 
IVo TBC (* = -l) 



[8, 7] [6, 5] [4,3] [2,1] = 
[7,6] [5,4] [3, 2] {1,8} 



+^o TBC (<5 



+1) 



-V>o TBC (5 = -l) 



(8) 
(9) 



we see that i[)q BC (5 = ±1) have different parity eigenvalues P = ±1. Then the states 
4)q BC (5 > 0) and ipj (5 < 0) have different parities P since they are adiabatically connected 
to ipo ($ = +1) an( A ipQ BC (S = +1), respectively. The Hamiltonian (2) is P-invariant, which 
means there is no mixing between the states with different P. Thus the lowest P = +1 state 
and P = — 1 state crosses with each other at the quantum phase transition point 5 = 0, as 
demonstrated in Fig. 3. We note that ipo(5 > 0) and tpo(5 < 0) have same eigenvalues P for the 
PBC case. 





Figure 2. Schematic behaviors of i/jq (S > 

0) and Vo BC (<5 < °) under the PBC - The y 
do not cross with each other (not broken 
lines). They smoothly changes with each 
other because the mixing of these two states 
occurs. 



Figure 3. Schematic behaviors of ip FBC (S > 

0) and Vo BC ( 5 < °) under the TBC - The y 
cross with each other at the transition point 
5 = 0. The mixing of these two states does not 
occur because they are protected by different 
P. 



3. S = 2 case 

Here we explain the S = 2 case, considering N = 4 system for simplicity. Under the PBC the 
Haldane state and the ID state corresponding to the valence bond pictures are, respectively 

^ BC = [1,2] 2 [2,3] 2 [3,4] 2 [4,1] 2 (10) 
^ P D BC = [1,2][2,3][3,4][4,1] (11) 



where two valence bonds between the ith and jth spins are abbreviated as [i,j] 2 - Although 
these wave functions are not exact ones in general, the exact wave functions for the Haldane 
and ID states are adiabatically connected to above two wave functions respectively. When we 
operate X on these states, we obtain 

J^ BC = [4,3] 2 [3,2] 2 [2,1] 2 [1,4] 2 =^ BC (12) 
X^ P D BC = [4, 3] [3, 2] [2, 1] [1,4]= Vi P D BC (13) 

Then, P = +1 for both the Haldane state and the ID state. Since there are no valence bonds 
in the LD state, it is clear P = +1 for the LD state. We cannot distinguish these three states 
under the PBC. 




Figure 4. Haldane state under the PBC. Figure 5. ID state under the PBC. 



Under the TBC the Haldane state and the ID state corresponding to the valence bond picture 
are, respectively 

^S BC = [1,2] 2 [2,3] 2 [3,4] 2 {4,1} 2 (14) 
^™ c = [1,2] [2, 3] [3, 4] {4,1} (15) 

When we operate X on these states, we obtain 

X^ BC = [4,3] 2 [3,2] 2 [2,1] 2 {1,4} 2 =^ BC (16) 
X^ T D BC = [4, 3] [3, 2] [2, 1] {1, 4} = -^ T D BC (17) 

Then, P = +1 for the Haldane state and P = — 1 for the ID state. Since there are no valence 
bonds in the LD state, it is clear P = +1 for the LD state. Thus we can distinguish the ID 
state from the Haldane and LD states by the eigenvalue P. Namely, under the TBC, if the 
lowest eigenstate is P = —1, the ground state is the ID state. On the other hand, if the lowest 
eigenstate is P = +1, the ground state is the Haldane state or the LD state. 

There may arise a question how to distinguish the Haldane state and the LD state. 
Recently Pollmann, Berg, Turner and Oshikawa [9] stated that the Haldane state is essentially 
indistinguishable from the LD state. They also constructed a one-parameter matrix product 
state which interpolates the Haldane state and the LD state without any quantum phase 
transition. Our results [1] also indicate that the Haldane state and the LD state belong to 
the same phase, as shown later in Figs. 10 and 11. 

4. Ground state phase diagram of an S = 2 chain 

From the discussion of the previous section, we should compare the two energy levels Ej BC {M = 
0,P = +1) and E^ BC (M = 0,P = -1), where E^ BC (M,P) is the lowest energy in the 
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subspace of magnetization M = J2j S'j and parity P under the TBC. To obtain the phase 
diagram, we further check whether the ground state is gapless (XY) state or gapped state 
(Haldane, ID or LD). Nomura and Kitazawa [7] showed that the ground state is gapless when 
E PBC (M = 2) < Ej BC (M = 0,P), where E PBC (M = 2) is the lowest energy in the subspace 
M = 2 under the PBC. This condition was obtained through the effective Hamiltonian and 
the renormalization group method, unfortunately for which we have no physical and intuitive 
explanation. Anyway, we have to compare three energy levels, Eq BC (M = 0,P = +1), 
Ej BC (M = 0,P = -1) and E PBC (M = 2). Namely, the ground state is the Haldane/LD (H/LD) 
state, the ID state or the XY state according as Pj BC (M = 0, P = +1), Ej BC (M = 0,P = -1) 
or E PBC (M = 2) is the lowest. 
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Figure 8. Behavior of AEq BC = Eq BC (M - 
0,P = +l)-^ rBC (M = 0,P = -l). 



PBC-TBC 



Figure 9. Behavior of AE 
E PBC (M = 2)-Ej BC (M = 0,P =+l). 



Figure 8 shows AE^ BC = P rBC (M = 0,P = +1) - P rBC (M = 0,P = -1) function 
of A when D = 1.82 obtained by the numerical diagonalization of 12 spin systems. From this 
figure we can determine the phase boundary between the Haldane/LD phase and the ID phase. 
We note that E PBC (M = 2) is higher than P rBC (M = 0, P = ±1) in this region. The behavior 
of AP BBC - TBC = E PBC (M = 2) - P rBC (M = 0,P = +1) is shown in Fig. 9, from which 
the phase boundary between the XY phase and Haldane/LD phases is determined. The energy 
P rBC (M = 0, P = -1) is higher in this region. Examples of the size dependence of the critical 
values of A are summarized in Table 1, from which we see that the finite-size effects are not 
so serious in our LS anayses. Since the transition between the Haldane/LD phase and the 
Neel phase is expected to of the Ising type, we use the phenomenological renormalization group 
method for the numerical analysis [1]. 

Our final phase diagram is shown in Figs. 10 and 11. The remarkable natures of the phase 



Table 1. Examples of critical values of A in the case of D = 1.82 obtained by the LS method. 
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2.17687 


2.14971 


2.08931 
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2.22262 


2.16106 


2.12246 


10 


2.23529 


2.16527 


2.13607 


12 


2.23793 


2.16639 


2.14265 


oo 


2.241 ±0.001 


2.167 ±0.001 


2.156 ±0.001 




Figure 10. Phase diagram of an S = 2 Figure 11. Enlarged phase diagram 
chain. H and N are the abbreviations of near the ID phase. 
Haldane and Neel, respectively. 

diagram are: (1) there exists the ID phase which was predicted by Oshikawa in 1992 and has 
been believed to be absent for a long time; (2) the Haldane state and the LD state belong to 
the same phase. 
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